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We report on experimental measurements of the growth of 
regular domains evolving from an irregular pattern in elec- 
troconvection. The late-time growth of the domains is con- 
sistent with the size of the domains scaling as t n . We use 
two isotropic measurements of the domain size: the structure 
factor and the domain wall length. Measurements using the 
structure factor are consistent with £ 1//5 growth. Measure- 
ments using the domain wall length are consistent with £ 1//4 
growth. One source of this discrepancy is the fact that the dis- 
tribution of local wavenumbers is approximately independent 
of the domain size. In addition, we measure the anisotropy of 
the growing domains. 
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There are many situations where a system experiences 
a rapid change of an external parameter, or quench, such 
that the state of the system after the quench is not an 
equilibrium or steady-state phase. Domains of the new 
equilibrium phase form, and the subsequent growth of 
these domains, or coarsening, is often characterized by a 
single scale size for the domains that follows a power-law 
growth. There has been a great deal of both theoreti- 
cal and experimental study of this process for systems 
where both the initial state prior to the quench and fi- 
nal state after the quench are thermodynamic equilib- 
rium states at finite temperature [Q. We are interested 
in the analogous process for systems that are driven out 
of equilibrium. For such systems, neither the initial nor 
the final state is in thermodynamic equilibrium; however, 
they are steady-states of the system. The first situation 
will be referred to as a thermodynamic system and the 
second situation as a driven system. For driven systems, 
studies of model equations suggest power-law scaling of 
the late-time domain growth; however, the value of the 
growth exponent depends on the measurement scheme 
[|]-[| . Experimentally, the coarsening of random patterns 
after a ramp in the control pattern has been studied, but 
growth exponents for domain size were not measured [[| . 
Thermodynamic systems with stripes have been studied 
experimentally using block copolymers ||. 

In spatially extended systems that are driven out of 
equilibrium, there is generally a transition from a uni- 
form state to a periodic, or "stripe", state at a criti- 



cal value R c of the external control parameter R [Q. A 
quench corresponds to a rapid change of R. For val- 
ues reasonably close to R c , these systems are often well 
described by model equations that are similar to, and 
in some cases identical to, the equations used to study 
coarsening in thermodynamic systems 0. However, in 
a driven system, there is generally no equivalent of a 
free-energy that governs the growth of the domains, and 
the periodic structure complicates the dynamics. There- 
fore, the question of how the ordering proceeds in driven 
systems, and the differences and similarities with ther- 
modynamic systems, is one of great interest. 

Simulations of potential ||,[| and nonpotential || 
forms of the Swift-Hohenberg equation have been used 
to study quenches in driven systems [||-f§. In both cases, 
characterization of the growth by the structure function 
S(k) suggests a length scale that grows as t 1 ^ HQ- In 
contrast, the growth exponent determined from the ori- 
entational correlation function is consistent with 1/4 [|]]4|] 
for potential dynamics and with 1/2 ^ for nonpotential 
dynamics. There is no explanation of the discrepancy 
between the S(k) and orientational correlation function 
measurements; however, the length scales determined by 
these measures do not have the same immediate interpre- 
tation The orientational correlation function results 
agree with experiments in block copolymers 

We have made experimental measurements of coarsen- 
ing in an anisotropic, driven system: electroconvection in 
a nematic liquid crystal [§ . A nematic liquid crystal is a 
fluid in which the molecules align on average along a par- 
ticular axis, referred to as the director ||, which we take 
as the x-axis. Because the system is anisotropic, there 
are only two possible orientations of the domains. This 
is fundamentally different than previous simulations and 
experiments. One expects a different scale size parallel 
and perpendicular to the preferred direction in the sys- 
tem. Our measurement of these scale sizes show that not 
only are the magnitudes of the length scales different, but 
they also coarsen with different growth exponents. We 
also measure isotropic properties of the domain growth 
and find that they are in surprising agreement with the 
predictions of the Swift-Hohenberg simulations ||,[| . Fi- 
nally, our measurements show that the domain size and 
wavenumber variation have different growth exponents, 
which explains the discrepancy between the S(k) and ori- 
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entational correlation function measurements. 

For electroconvection, the liquid crystal is placed be- 
tween two glass plates with the director aligned parallel 
to the plates and along a single axis. The liquid crystal 
is doped with ionic impurities. An ac voltage is applied 
perpendicular to the plates. Above a critical value of the 
applied voltage V c , there is a transition to a state that 
consists of convection rolls with a corresponding periodic 
variation of the director and charge density. We studied 
oblique rolls, i.e. patterns that have a nonzero angle 
between the wavevector and the alignment of the undis- 
torted director. Oblique rolls with the same wave number 
k but at (zig rolls) and —6 (zag rolls) are degenerate. 
The initial transition in this system is to a pattern that 
consists of the superposition of four modes: right- and 
left-traveling zig and zag rolls. The interaction of these 
four modes leads to irregular spatial and temporal varia- 
tions of the amplitudes of each of the modes [flOfl , i.e. spa- 
tiotemporal chaos For traveling rolls, a sufficiently 

large modulation of the amplitude of the applied voltage 
at twice the intrinsic frequency of the pattern stabilizes 
standing rolls |]l^|-|l4[. For our system, either standing 
zig or standing zag rolls are stabilized, and the stabilized 
pattern exhibits regular temporal dynamics |l5|] . In our 
experiments, a quench corresponded to a rapid change of 
the modulation amplitude at a fixed value of V rms . As 
standing waves can be stabilized both below and above 
V C1 two types of quenches are possible. Below V c , the 
dynamics of the standing waves are potential and 
above V c , they are nonpotential [ fl7| |. Therefore, based on 
the results of Ref. || , we expected two different growth 
exponents. Also, below V C1 domains of standing zig and 
zag rolls must first form before the system coarsens, as 
the initial state is spatially uniform. This process is illus- 
trated in Fig. 1(a) - (d). Above V c , domains of zig and 
zag rolls are already present after the quench, and the 
coarsening proceeds from this initial distribution. This 
process is illustrated in Fig. 1(e) - (h). 

The details of the experimental apparatus are de- 
scribed in Ref. (l5| . Commercial cells |l8| with a thick- 
ness of 23 fim and 1 cm 2 electrodes were used. The 
sample temperature was maintened at 50.0 ± 0.002 °C. 
The patterns were observed from above using a modified 
shadowgraph setup [[l9| that effectively eliminated the 
well-known nonlinear effects in the shadowgraph image 
[|o[. An ac voltage, V(t) = [V Q + V m sm(u m t)] sin(cc^), 
was applied across the sample, with Ud/^n = 25 Hz. 
There are two relevant dimensionless control parameters: 
e = (Vo/Vc) 2 — 1 and b = V m /V c . Here V c is the critical 
voltage for the onset of a pattern when Vm = 0. There 
was a small, linear drift in V c \ therefore, before and after 
each quench, V c and ojh were determined. Typical values 
were V c = 15 V rms and Uh/Zn = 0.5 Hz. 

The system was equilibrated at b = and either e = 
—0.03 or e = 0.03 for five minutes. For e = —0.03 (e = 
0.03), the transition to standing waves occurs at b = 0.03 



(b = 0.013) Jl^| . For e = —0.03, we used a jump from 
b = to b = 0.05, and for e = 0.03, we used a jump from 
b = to b = 0.04. In both cases, uo m = 2^, where uoh is 
the frequency of the pattern at onset, the Hopf frequency. 
A third quench was done at e = 0.03 with co m = 2.\uj\ 1 . 
Immediately after a jump in 6, a series of 128 images 
was obtained. An image was taken every 5 th cycle of the 
modulation, or roughly once every 5 seconds, triggered by 
the applied voltage to occur at T/4, where T is the period 
of the modulation. In order to resolve the individual rolls, 
a 1.35 mm by 1.35 mm section of the sample was imaged. 
At the end of a typical time series, the domain size was 
on the order of our observation window; however, given 
the size of the entire sample, we were not limited by 
finite size effects. Unless otherwise noted, time is scaled 
by the director relaxation time, which for our system is 
0.2 s. For e = 0.03, we performed 20 quenches, and for 
e = —0.03, we performed 18 quenches. The domain wall 
length, the local wavenumber distribution, Aq x , 

and Aq y were computed for each image in the time series 
for a given quench. (Aq x and Aq y are defined below.) 
The values at each time step were averaged over all the 
quenches for a given e. These averaged values were used 
to compute the growth exponents, which did not change 
significantly after averaging 10 quenches. 

If a single scale length is sufficient to describe the do- 
main growth, both the total domain wall length and the 
width of the power spectrum would exhibit the same 
growth exponent. Because, the structure factor con- 
tains information about the range of local wavenumbers 
within a domain, we measured the local wavenumber 
using the method described in Ref. pi] . Briefly, this 
method involves calculating the x-component of k from 
\k x \ 2 = — [d 2 u(x.)]/u(x), with a similar calculation for the 
y-component. Here u{x) is the image of the pattern. We 
used the central 0.82 mm by 0.82 mm region of the image 
in the determination of the spread Ak of the magnitude 
k of the local wavevector k. We defined Ak as the square 
root of the second moment of the distribution of k about 

To determine the normalized domain- wall length L, 
the local orientation was determined from the sign of 
[d X yu(x.)]/u(x). Regions with a negative (positive) value 
of this ratio corresponded to zig (zag) rolls and were as- 
signed a grayscale of (255). The resulting image was 
smoothed, producing an image in which domain walls 
had a value close to 128. Applying a threshold produced 
images in which zig regions had a value of 0, domain walls 
a value of 128, and zag rolls a value of 255. The domain 
wall "length" was taken as the total number of pixels of 
value 128 normalized by the total number of pixels in the 
image. Images from a typical sequence, after processing, 
are shown in Fig. 1. The processed images were also used 
to study the anisotropy of the growth. Since the wave- 
length has been factored out of the processed images, the 
width of the central peak in the power spectrum 5(q) of 
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one of these images corresponds to the length scale of the 
domains. Defining Aq x (Aq y ) as the second moment of 
Qx (Qy) about q = of *S(q), with q x and q y the x- and 
y-components of q, provides a measure of the inverse of 
the correlation length in the x and y directions. 
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FIG. 1. Processed images illustrating the domain growth 
for a quench at e = —0.03 (images (a) - (d)) and for a quench 
at e = 0.03 (images (e) - (h)). The images cover an area of 
1.35 mm by 1.35 mm, and the processing method is described 
in detail in the text. The black regions are areas of zig rolls, 
and the white areas are regions of zag rolls. Image (a) / (f) was 
taken 5 s after the quench, and the subsequent images are all 
160 s apart. 

The measurements based on S(k) used the same 
method as described in Ref. ||. The structure factor 
(the square of the Fourier transform) was averaged over 
all angles. The relevant peak in the resulting S(k) was fit 
to a Lorentzian squared, and the width Sk was defined as 
the half width at half height. This method is an isotropic 
measure of the domain growth and is used to ensure that 
our results are directly comparable to Ref [|| . 



i — | — i — | — i — | — i — | — i — | — I I — i — | — i — | — i — | — i — | — i — r 




1.0 1.5 2.0 2.5 3.0 3.5 1.0 1.5 2.0 2.5 3.0 3.5 

log (t) log (t) 

FIG. 2. This shows a plot of log(L) versus log(t). Here L 
is the total length of the domain wall in the region of study. 
Plot (a) is for the quench at e = —0.03, and plot (b) is for the 
quench at e = 0.03. The solid lines are linear fits to the data. 
The line in (a) has a slope of —0.240 =b 0.004, and the line in 
(b) has a slope of —0.243 =b 0.003. In (b), the open symbols 
are the results for Aq x (circles) and Aq y (squares). 

Figures 2a and 2b show the results for the domain wall 
length L for the e = —0.03 and e = 0.03 quenches, re- 
spectively. Here log(L) is plotted versus log(t). For the 
quench at e = 0.03, the behavior is consistent with power 
law scaling essentially immediately after the quench. For 



e = —0.03, the system is not consistent with power law 
scaling until log(t) « 2.5. This difference is reasonable 
given that the domains must first form for the quench 
at e = —0.03. Also, the scaling occurs in both sys- 
tems at roughly the same scale size for the domains, 
log(L) w —1.2. The decay of the domain-wall length is 
consistent with scaling as £ -1 / 4 . The solid line in Fig. 2a 
is a fit over the range 2.5 < log(t) < 3.5 and has a slope 
of —0.24. The solid line in Fig. 2b is a fit over the range 
2.0 < log(t) < 3.2 and also has a slope of —0.24. 

Figure 2b also shows the results for Aq x (open circles) 
and Aq y (open squares). As is suggested by the images 
in Fig. 1, the domains tend to be larger along the direc- 
tor (x-direction). This is confirmed by the relative mag- 
nitudes of Aq x and Aq y . Also, the two lengths coarsen 
with different exponents, where the growth exponent per- 
pendicular to the director is consistent with the exponent 
measured using the domain-wall length. A similar result 
holds for the e = —0.03 quench, with the same delay in 
the onset of scaling as seen in the wall length. 
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FIG. 3. This shows a plot of \og(Sk) versus log(t) using the 
left hand axis and solid symbols, and a plot of log(Afc) versus 
log(t) using the right hand axis and open symbols. Plot (a) 
is for the quench at e = —0.03, and plot (b) is for the quench 
at e = 0.03. The solid lines are linear fits to the solid data. 
The line in (a) has a slope of —0.170 =b 0.003, and the line in 
(b) has a slope of -0.165 ± 0.002. 

Figures 3a and 3b show \og(Sk) and log(Afc) versus 
log(t) for the quenches at e = —0.03 and e = 0.03, respec- 
tively. The behavior is consistent with power-law scaling 
almost immediately for e = 0.03 and at later times for 
e = —0.03. The scaling sets in for both quenches at ap- 
proximately the same domain size. For both quenches, 
the scaling of Sk(t) is consistent with £ -1 / 5 decay, as 
found in simulations HQ. For e = —0.03, the solid line is 
a fit over the region 2.5 < log(t) < 3.5 and has a slope of 
-0.17. For e = 0.03, the solid line is a fit over the region 
2.0 < log(t) < 3.2 and has a slope of -0.16. However, for 
both cases, the variation in the local wavenumber Ak is 



3 



a significant fraction of 5k in the possible scaling regime 
and is roughly constant in time. Therefore, 5k is not an 
accurate measure of the growth of the domain size, but a 
complicated convolution of the domain size and Ak. Ac- 
counting for this effect, measurements of the coarsening 
parallel and perpendicular to the director based on indi- 
vidual peaks of S(k) are consistent with the anisotropy 
determined from *5(q) of the processed images. 

Despite the different symmetries, the isotropic mea- 
sures of the domain growth in electroconvection agree 
well with the simulations of the Swift-Hohenberg model 
and experiments in block copolymers ||. Under- 
standing this agreement will require further work. In 
particular, the work with block copolymers attributes the 
exponent of 1/4 to the dynamics of topological defects 
||. Because of the anisotropy, the relevant defect dy- 
namics in electroconvection are different. For example, 
there are no regions where the roll orientation changes 
continuously, only sharp domain walls. 



1.08 
1.06 
1.04 

o 

1.02 

> 

1.00 
0.98 



_o 

.o 



_l_ 



_l_ 



100 200 

time (s) 



300 



FIG. 4. This shows a plot of k a ve/k versus time. The 
triangles are for the quench at e = 0.03 with cj m = 2.0cJh,, and 
the circles are for the quench at e = 0.03 with co m = 2.1uoh- 
The lines are a guide to the eye. 

One discrepancy is our observation of a growth ex- 
ponent of 1/4 for both e < and e > 0. For e > 0, 
we expected an exponent of 1/2 [||. One possible ex- 
planation of the discrepancy is the potential existence of 
long-time transients. An analysis of the Swift-Hohenberg 
equation for patterns with k = k Q suggests a long-time 
transient regime for which the growth exponent is 1/4 
• In our system, the dispersion relation fixes the av- 
erage wavenumber k ave = k Q for uo m = 2coh- To vary 
kave, we performed a quench at co m = 2.1^, for which 
k ave = 1.07/c o . The quench was at e = 0.03 to a value 
of b = 0.04. Measurements of L were consistent with a 
growth exponent of 1/4. However, the average wavenum- 
ber exhibited a slow approach to a steady state (see Fig. 
4), suggesting that this growth exponent also corresponds 
to a transient regime. Because there is no significant evo- 
lution of the wavenumber spread, we plan to study larger 
regions of the sample, without resolving the individual 
rolls, for longer times by exploiting optical properties of 



the patterns. This will allow for the possibility of observ- 
ing a crossover to a different growth exponent. 

The authors acknowledge funding from NSF grant 
DMR-9975479. L. Purvis was supported by NSF 
Research Experience for Undergraduates grant PHY- 
9988066. M. Dennin would like to thank Michael Cross 
and Herman Riecke for useful discussions. 



[1] A. J. Bray, Advances in Physics 43, 357 (1994). 

[2] K. R. Elder, J. Vinals, and M. Grant, Phys. Rev. Lett. 

68, 3024 (1992). 
[3] M. C. Cross and D. I. Meiron, Phys. Rev. Lett. 75, 2152 

(1995). 

[4] J. J. Christensen and A. J. Bray, Phys. Rev. E 58, 5364 

(1998) . 

[5] C. W. Meyer, G. Ahlers, D. S. Cannell, Phys. Rev. A 44, 
2514 (1991). 

[6] C. Harrison, D. H. Adamson, Z. Cheng, J. Sebastian, S. 

Sethuraman, D. Huse, R. A. Register, and P. M. Chaikin, 

Science 290, 1558 (2000). 
[7] For reviews of pattern formation, see M. C. Cross and P. 

C. Hohenberg, Rev. Mod. Phys. 65, 851 (1993), and J. 

P. Gollub and J. S. Langer, Rev. Mod. Phys. 71, s396 

(1999) . 

[8] Review articles on electroconvection can be found in I. 
Rehberg, B. L. Winkler, M. de la Torre Juarez, S. Rase- 
nat, and W. Schopf, Festkorperprobleme 29, 35 (1989); 
S. Kai and W. Zimmermann, Prog. Theor. Phys. Suppl. 
99, 458 (1989); and L. Kramer and W. Pesch, Annu. Rev. 
Fluid Mec. 27, 515 (1995). 

[9] P. G. de Gennes, The Physics of Liquid Crystals (Clare- 
don Press, Oxford, 1974); S. Chandrasekhar, Liquid 
Crystals (Cambridge University Press, Cambridge, Eng- 
land, 1992). 

[10] M. Dennin, G. Ahlers, and D. S. Cannell, Science 272, 
388 (1996). 

[11] For a recent review of experimental examples of spa- 
tiotemporal chaos, see G. Ahlers, Physica A 249, 18 
(1998). 

[12] H. Riecke, J. D. Crawford, and E. Knobloch, Phys. Rev. 

Lett. 61, 1942 (1988). 
[13] D. Walgraef, Europhys. Lett. 7, 485 (1988). 
[14] I. Rehberg, S. Rasenat, J. Fineberg, M. de la Torre 

Juaurez, and V. Steinberg, Phys. Rev. Lett. 61, 2449 

(1988). 

[15] M. Dennin, Phys. Rev. E, (2000). 

[16] H. Riecke, Europhys. Lett. 11, 213 (1990). 

[17] H. Riecke, private communication. 

[18] E.H.C. CO., Ltd., 1164 Hino, Hino-shi, Tokyo, Japan. 

[19] H. Amm, R. Stannarius, and A. G. Rossberg, Physica D 

126, 171 (1999). 
[20] S. Rasenat, G. Hartung, B. L. Winkler, and I. Rehberg, 

Experiments in Fluids 7, 412 (1989). 
[21] D. A. Egolf, I. V. Melnikov, and E. Bodenschatz, Phys. 

Rev. Lett. 80, 3228 (1998). 



4 



